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A  TORQUE  MAGNETOMETER  FOR  CRYSTALLINE 
ANISOTROPY  MEASUREMENTS  * 
by 

Robert  P.  Hunt 

Laboratory  for  Insulation  Research 
Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 

Abstract:  An  instrument  designed  to  measure  the  torque  of  magnetic  samples 
automatically  and  continuously  as  a  function  of  the  angle  of  an  applied 
magnetic  field  with  respect  to  the  sample  is  described  in  detail.  Gen¬ 
eral  operation  of  the  equipment  is  outlined,  and  the  system  is  analyzed 
on  the  basis  of  feedback  theory.  Results  of  a  temperature  run  on  mag¬ 
netite  and  yttrium  iron  garnet  are  shown,  along  with  a  room-temperature 
determination  of  the  anisotropy  of  nickel.  Measurements  have  also  been 
made  on  the  induced  anisotropy  of  polycrystalline  yttrium  iron  garnet 
and  associated  relaxation  phenomena. 

Magnetocrystalline  anisotropy  energy  will  generally  be  an  important  con¬ 
sideration  in  any  detailed  study  of  magnetic  materials.^  This  paper  describes 
an  instrument  which  facilitates  the  rapid  determination  of  the  anisotropy  energy 
of  magnetic  single  crystals. 

There  are  three  basic  means  by  which  the  anisotropy  may  be  experi¬ 
mentally  determined.  The  first  is  an  energy  evaluation  from  the  magnetization 

curves  taken  with  the  field  at  various  directions  to  the  crystalline  axes.  This 
* 

Based  on  a  thesis  submitted  in  partial  fulfillment  of  the  requirements  for  the 
Degree  of  Master  of  Science  in  Electrical  Engineering  at  the  Massachusetts 
Institute  of  Technology, 

1)  R.  M.  Bozorth,  "Ferromagnetism,  "  D.  Van  Nostrand  Co.  ,  New  York, 

1955. 
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is  a  lengthy  and  tedious  procedure.  The  second  is  based  on  ferromagnetic 
resonance  experiments.  This  method  cannot  be  used  on  materials  exhibiting 
a  skin  effect;  also  the  anisotropy  of  some  materials  may  have  an  appreciable 
frequency  dependence. 

The  final  method  makes  use  of  the  directional  dependence  of  the  anisotropy 
energy  to  generate  a  torque  which  can  be  measured  by  a  torsion  balance.  Here 
the  anisotropy  is  determined  under  static  conditions,  but  again  the  procedure  is 
lengthy.  To  overcome  this  difficulty  we  have  designed  a  torque  magnetometer 
that  measures  the  crystalline  torque  automatically  and  continuously  as  a  function 
of  the  direction  of  the  applied  field. 

The  instrument  is  capable  of  measuring  torques  from  5  to  20,  000  dyne -cm 
over  several  ranges  with  a  full-scale  accuracy  of  2%.  The  torque  is  plotted  on 
the  Y  axis  of  an  X-Y  recorder.  The  X  axis  plots  the  angle  of  the  applied  field, 
which  may  be  rotated  about  the  crystal.  Low  torque  measurements  are  some¬ 
what  hampered  by  background  vibrations.  An  entire  torque  curve  may  be  taken  in 
about  one  minute.  Equipment  for  controlling  temperatures  from  liquid-nitrogen 
to  400°C  is  provided. 


The  Torque  Magnetometer 
Description  of  Equipment 

The  instrument  is  basically  a  self-balancing  torsion  balance  which  sub¬ 
stitutes  an  electronic  feedback  system  for  the  torsion  fiber.  The  essentials  of 
the  unit  are  shown  in  Fig.  1. 

The  field  applied  to  the  crystal  is  derived  from  a  Varian  six-inch  magnet, 
mounted  on  a  gun-sight  mount,  so  that  the  field  may  be  rotated  about  a  central  axis 
at  the  midpoint  of  the  magnet's  gap.  A  ten-turn  potentiometer  ganged  to  the  mag¬ 
net  provides  the  electrical  indication  of  the  magnet's  position,  which  drives  the 


2)  L.  R.  Bickford,  Phys.  Rev.  76,  137  (1949);  78,  449  (1950). 


-  3  - 


Permanent 

magnet 


Light  source 


Balancing  y  ft*-' 
coil  y  '' 

Mirror  ✓ 

Teflon  bushing 


Balancing 

amplifier 


— o 

<  To  recorder 

<  y  axis 

1 - o 


Fig.  1.  Torque  magnetometer . 


X  axis  of  an  X-Y  recorder. 

The  sample  is  placed  at  the  midpoint  of  the  air  gap  and  is  suspended  by 

a  shaft,  cn  the  upper  end  of  which  is  attached  a  d'Arsonval  type  coil  assembly, 

3) 

similar  to  that  described  by  Penoyer,  '  which  is  situated  between  the  poles  of  a 
small  magnetron  magnet  (Fig.  2).  The  coil  and  suspension  mechanism  pivot  on 
a  jewel  bearing;  the  shaft  is  guided  by  a  Teflon  knife-edge  bushing  to  reduce 
friction. 

Below  the  coil  assembly  a  small  galvanometer  mirror  is  fixed  to  the  shaft. 
A  stationary  light  source  focuses  its  beam  on  the  mirror  and  is  reflected 
onto  the  cathodes  of  a  dual  phototube,  which  is  incorporated  in  the  balancing 

3)  R.  F.  Penoyer,  "An  Automatic  Torque  Balance,  "  IBM  Technical  Publication, 
IBMCorp.,  Poughkeepsie,  New  York,  1957. 
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amplifier.  The  beam  is  so  adjusted  that  normally  an  equal  portion  of  light  flux 
falls  on  each  cathode. 

When  a  torque  is  applied  to  the  sample  shaft,  the  beam  is  displaced  to 
one  side  of  the  phototube.  The  latter  "senses"  the  displacement  of  the  shaft  and 
produces  an  error  voltage  proportional  in  magnitude  and  sign  to  the  displace¬ 
ment.  The  amplifier  converts  the  error  signal  to  a  d-c  restoring  current  which 
is  fed  to  the  balancing  coil.  The  torque  developed  by  the  balancing  coil  is  equal 
but  opposite  to  the  applied  torque  of  the  crystal.  A  direct  measure  of  this  ap¬ 
plied  torque  is  the  current,  which  is  then  plotted  on  the  Y  axis  of  the  X-Y  re¬ 
corder.  The  entire  system  is  a  closed-loop,  self-balancing,  feedback-controlled 
unit. 

With  this  apparatus  the  torque  exerted  by  a  crystal  sample  versus  the 
angle  of  the  applied  field  can  be  plotted  easily  and  quickly  (cf.  Fig.  19,  p.  28). 
The  angle  of  the  field  may  be  readily  related  to  the  crystallographic  axes,  which 
are  assumed  to  be  known  a  priori. 

Design  Criteria 

Anisotropy  measurements  are  generally  made  on  single -crystal  disks  of 

about  1  cm  or  less  or  on  small  spherical  samples  of  a  few  millimeters  diameter, 

cemented  to  the  end  of  the  sample  shaft.  Care  must  be  taken  to  avoid  attaching 
*  _ 

Orientation  is  determined  by  X  rays  prior  to  experiment. 
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the  sample  so  rigidly  (as  with  a  spring-loaded  clamp)  that  the  crystal  is  under 

4) 

a  strain  since  this  would  affect  the  anisotropy  behavior.  '  The  field  gap  of  the 
Varian  magnet  is  two  inches,  allowing  plenty  of  room  for  a  Dewar  or  heating 
furnace. 

A  survey  of  the  materials  that  had  been  previously  measured  indicated 

x 

that  the  torque  magnetometer  should  be  capable  of  measuring  torques  from  a 
4 

few  up  to  about  10  dyne -cm.  Together  with  a  knowledge  of  the  amplifier - 
current  output  and  the  magnetron  magnet's  field  B^,  this  provides  the  infor¬ 
mation  necessary  for  the  design  of  the  balance  coil.  The  Varian  magnet  must 
be  capable  of  saturating  the  crystal  in  question  to  avoid  analytical  compli¬ 
cations. 

Since  the  device  is  a  feedback  system,  the  frequency  response  must  be 
considered.  It  clearly  must  be  flat  over  the  region  of  interest.  This  region 
may  be  determined  from  the  field-angie  dependence  of  the  torque  and  the 
angular  velocity  of  rotation  of  the  magnet. 

Bozorth^  showed  that  the  most  "rapid"  significant  variation  of  torque 
for  the  cubic  and  hexagonal  structure  has  a  60  dependence.  Hence  we  will 
assume  this  to  be  our  highest-order  effect.  0  is  determined  by  the  angular 
velocity  of  the  field,  w  ,  i.  e.  , 


0 


=  CO 

m 


t, 


(1) 


where  ^  2it/15  radians /sec  and  is  limited  by  the  magnet's  drive  motor. 
Hence,  the  maximum  radian  frequency  to  be  encountered  is  6com  or  12ir/l5 
radians/sec,  and  the  unit  must  be  linear  up  to  this  value. 


4)  D.  Gugan  and  G.  Rowlands,  Proc.  Phys.  Soc.  (London)  72,  207  (1958). 

5)  R.  M.  Bozorth,  Phys.  Rev.  50,  1076  (1936). 
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System  Analysis 
Block-Diagram  Representation 

From  Fig.  1  and  the  description  of  the  equipment,  a  general  block  dia¬ 
gram  (Fig.  3)  may  be  formed  for  the  system.  For  convenience  we  shall  work 


Fig.  3.  General  block  diagram. 

mainly  in  the  frequency  domain  making  use  of  the  Laplace  transform  theory.  ^ 
L.(s)  is  defined  as  the  Laplace  transform  of  the  input  torque  function;  L  (s)  is 
an  error  torque  given  by  the  difference  between  input  torque  L^(s)  and  output 
torque  Lq(s).  This  error  torque  acts  on  the  shaft  to  produce  a  small  angular 
deflection  4>(s)  which  generates,  through  the  optical  system  and  phototube  sen¬ 
sor,  an  error  voltage  E(s).  E(s)  is  converted  by  the  amplifier  to  a  current  l(s), 
which  in  turn  produces  the  output  torque  Lq(s)  through  the  balance  coil.  In  each 
case  the  parameters  may  be  inversely  transformed  to  find  the  corresponding 
time -domain  function;  i.  e.  ,  L^(a) -<— ►  L^(t),  etc. 

We  now  must  formulate  the  transfer  function  for  the  various  blocks. 

Mechanical  Equation 

The  only  mechanical  equation  is  that  which  governs  the  shaft  motion.  If 


6)  J.  A.  Aseltine,  "Transformation  Methods  in  Linear  System  Analysis,  " 
McGraw-Hill  Book  Co.,  New.  York,  1958. 
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J  is  the  moment  of  inertia*  of  the  shaft,  a  the  damping  factor,  and  K  the  spring 
constant,  and  if  4>  is  the  angle  of  the  balance  coil's  area  vector  with  the  faces 
of  the  magnetron  magnet's  pole  pieces  (Fig.  2),  we  may  write  Newton's  law  for 


the  shaft  as 


T  d>  d<J>  . 

d7  31 


K+  =  L-.(t)  +  L£  -  L0(t). 


The  frictional  torque  was  found  to  be  10  newton-meters  (1  dyne  -cm)  and 
may  be  neglected.  Thus, 


_  d  <(>  dd> 

dtZ  ^ 


+  K4>  =  L.  -  L  . 
at  ^  i  o 


Transforming  to  the  frequency  domain  and  assuming  that  for  initial  conditions 
4>(0)  =  0  and  4>'(0)  =  0,  we  get 

Js2<|>(s)  +  as<|>(s)  +  K<j>(s)  =  L^(s)  -  I_o(s).  (4) 

Solution  for  <J>(s)  gives  a  relation  between  error  torque  and  angular  deflection: 


♦(»)  = 


L..(s)  -  L  (s) 
'  or  ' 


Le(s) 


Js  +  as  +  K  Js  +  as  +  K 


The  shaft  is  a  long,  hollow  brass  rod,  having  a  moment  of  inertia  given 


J  =  m(R2  +  R2)  ,  (6) 

where  m  is  the  mass,  and  and  the  inner  and  outer  radii,  respectively. 

The  damping  factor  a  arises  partly  from  friction  and  windage  which  can¬ 
not  be  calculated,  but  principally  from  the  action  of  the  balance -coil  form, 

which  is  a  rectangular  brass  loop  and  acts  as  a  shorted  turn  in  the  field  B  of 

m 

the  magnetron  magnet.  The  spring  constant  K  comes  primarily  from  the  tor¬ 
sion  of  the  wire  leads  connected  to  the  balance  coil  and  is  small  but  not  negligible. 


In  the  system  analysis  we  shall  use  mks  rationalized  units.  Some  of  the  cgs 
units  will  be  indicated  for  reference. 
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Optical  Gain  Factor 

The  optical  gain  factor  is  defined  as  the  output  voltage  of  the  sensing 
unit  per  radian  of  angle  rotated  by  the  shaft: 


optical  gain  =  ~ 


=  constant. 


(?) 


A  linear  variation  of  E(t)  with  <i>  is  based  on  the  assumption  that  the  output 
current  of  a  phototube  is  an  approximately  linear  function  of  the  amount  of  flux 
falling  on  its  cathode.  As  <$>  varies  linearly,  the  light  flux  on  one  section  of  the 
dual  phototube  increases  while  that  on  the  other  decreases  linearly.  This 
assumption  was  verified  by  experiment.  Hence, 


E(t)  =  Kj4>(t), 
E(s)  =  K^s). 


E{t)  is  physically  the  output  of  the  difference -amplifier  circuit  (Fig.  4). 


Phototube  circuit 


Fig.  4.  Circuit  diagram. 
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The  dual  phototube  forms  a  bridge  circuit  along  with  the  variable  1-megohm 
resistor  R^.  A  3-volt  400-cycle  signal  is  applied  across  the  phototube. 

When  <j>  =  0,  an  equal  light  flux  falls  on  each  section,  and  each  tube  has  the 
same  current.  The  signals  at  both  grids  of  the  12AX7  difference  amplifier  are 
equal,  and  the  difference -amplifier  output  is  zero.  As  4>  deviates  from  0,  a 
misbalance  in  the  bridge  occurs,  and  the  output  voltage  E(t)  is  a  measure  of  the 
difference  between  the  grid  signals  and  hence  of  4>.  E(t)  is  a  400-cycle  signal 
whose  magnitude  is  proportional  to  <}>.  The  phase  of  E(t)  (or  its  sign)  is 
0  degrees  for  <J>  >  0  and  180  degrees  for  <j>  <  0.  Another  approach  is  to  think  of 
the  phototube  circuit  as  a  chopper  system  acting  at  a  400 -cycle  rate. 

The  steps  used  in  arriving  at  Eqs.  8  are  shown  diagrammatically  in 
Fig.  5.  In  Fig.  5a,  E(t)  is  related  to  <J>  through  the  various  constants  of  the 
phototube  and  the  difference -amplifier  gain.  The  constants  may  be  lumped  into 
a  single  optical-gain  parameter  in  the  expression  for  E(t).  In  Fig.  5b,  the 
the  sine  multiplier  is  dropped,  since  it  is  effectively  canceled  out  in  the  de¬ 
modulator  circuit,  as  will  be  shown. 

Amplifier 

Since  E(t)  is  a  400-cycle  signal  modulated  with  the  essential  information 
about  <}>,  it  may  be  amplified  by  conventional  a-c  means.  This  is  the  function 
of  the  second  12AX7  in  Fig.  4.  The  signal  is  then  passed  through  a  transformer 
to  a  phase -sensitive  detector  where  the  400-cycle  signal  is  demodulated. 

After  detection  it  is  transmitted  through  a  filter  to  the  grids  of  a  6AS7  acting 
as  a  current  amplifier;  the  balance  coil  is  connected  across  the  cathodes  of  this 
tube.  The  6AG7  is  merely  an  oscillator  that  provides  the  400 -cycle  chopping 
and  detecting  voltage. 

* 


The  phototubes  rectify  the  current,  giving  half-wave  pulses.  This  does  not 
affect  the  above  analysis. 
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ea  =  ci  ■f"[<£c+<£(t)] 


ec(t) 


ebaCi4*[+c"^t)] 

c,  =  Phototube  constant 
g  =  Gain  of  difference  amplifier 
A  =  Area  of  the  two  phototube  cathodes 
Kt  =  Optical  gain  factor 
cu  =  Chopper  frequency 

(a) 


ec(t)  =c,gA<£(t)sina>t 
=  E(t)  sin  uJt 

E(t)  =  c,g<£(t)A 
=  K,<£(t) 


E(t) 

or  E(s)  =  K,<£(s) 
E(s) 


Fig.  5.  Phototube  circuit  operation:  (a)  functional 
diagram,  (b)  equivalent  representation. 


For  a  static  deflection  4>q  of  the  shaft,  the  output  of  the  phase -sensitive 

detector  is  a  d-c  voltage  whose  magnitude  is  a  measure  of  4>q.  If  the  voltage 

is  greater  than  0,  4»q  >  0,  and  vice  versa. 

In  the  frequency  domain  a  transfer -conductance  function  F(s)  for  the 

amplifier -demodulator -current  amplifier  system  is  defined  as 

Tp/.x  -  l(s)  _  output  current  ,  . 

E(  s)  input  voltage  ' 

The  a-c  amplifier,  since  it  operates  at  400  cycles,  provides  no  other  con¬ 
tribution  to  F(s)  than  a  straight  frequency-independent  gain  multiplier. 


Oscillator 


(a) 


(b) 


I(t)  =  g,g2  E(t)®a(t) 
=  K2E(t)®a(t) 


E(t) 


I(s) 

E(s) 


=  K2A(s) 

K2sts2 
"  (s+s,)(s+s2) 


Fig.  6.  Amplifier  detector: 


(a)  operation,  (b)  reduced  diagram. 


Similarly,  there  are  no  reactive  elements  in  the  d-c  current  amplifier.  Hence, 
any  poles  or  zeros  for  the  unit  must  originate  in  the  detector  circuit. 

The  operation  of  the  amplifier  is  shown  schematically  in  Fig.  6.  In 
Fig.  6b,  the  l/sin  tot  multiplier  has  been  omitted  since  it  is  effectively  canceled 
by  the  corresponding  multiplier  in  the  phototube  circuit.  The  various  tube  gain 
factors  have  been  concentrated  in  the  single  gain  factor  K^.  The  symbol  a(t) 
represents  the  impulse  response  of  the  filter  and  detector  circuit,  and  l(t)  is 
related  to  E(t)  in  the  time  domain  by  the  convolution  of  K2a(t)  with  E(t),  i.  e.  , 

I(t)  =  K2a(t)  ®  E(t)  =  K2  T00  a(T  )E(t  -  t  )d*  .  (9) 

Jo 

The  frequency-domain  equivalence  of  this  is  given  in  Fig.  6b,  where 

a(t)  ■«  ■>  A(s),  and  A(s)  is  determined  by  the  filter  and  detector  circuits.  It  is 

noted  that 

I(t)  I(s)  =  K2^>/’°°a(T)E(t  -T  )dx  =K2A(s)E(s).  (10) 

Jo 
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Fig.  7.  Detector -circuit  operation:  (a)  complete  cycle, 
(b)  first  half-cycle,  (c)  second  half-cycle. 


The  principal  elements  of  the  detector  circuit  are:  the  coupling  trans¬ 
former,  the  diodes,  and  the  filter  circuit.  The  diodes  switch  the  filter  from 
one  section  of  the  transformer  secondary  to  the  other  on  alternate  half-cycles. 
This  is  clarified  in  Fig.  7,  where  E  ^  is  a  large  400-cycle  reference  signal 
that  sweeps  the  diodes  of  each  ring  alternately  open  and  closed  on  each  half¬ 
cycle.  To  the  filter  circuit  the  diodes  appear  as  a  small  resistance  (effective¬ 
ly  the  forward  resistance  of  a  single  diode  in  series  with  the  transformer 
secondary).  For  analytical  purposes  the  circuit  may  be  simplified  to  that  of 
Fig.  8,  where  E  is  the  Thevenin  equivalent  voltage  output  of  the  transformer, 
rt  the  secondary  resistance  (one  section),  L.^  the  secondary  inductance,  r^ 
the  diode  forward  resistance,  R  an  exterior  resistor  (4400  ohms),  C  the  filter 
capacitance,  and  the  grid  resistor  of  the  6AS7.  Actually  r^  is  so  small  and 

R  so  large  that  both  may  be  dropped  from  consideration.  The  circuit  of  Fig.  9 

8 

represents  an  adequate  approximation  to  the  detector -filter  circuit.  This  is 
a  simple  R-L.-C  series  circuit  which  has  the  transfer  function 


eQ(s) 

w 


+  rt)s  r 


s2  + 


ctc’ 


'n\ 


where  c  is  a  constant  multiplier.  The  roots  of  the  denominator  of  Eq.  11  turn 
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where  F(s)  becomes  a  single  constant  gain  for  s  =  0  (i.e.,  is  defined 
as  the  static  gain). 
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Balance -Coil  Constant 

The  final  stage  of  the  system  is  the  conversion  of  the  output  current  of 
the  amplifier  to  a  torque.  This  is  accomplished  by  the  balance  coil  located 
between  the  poles  of  the  magnet'- on  magnet.  The  torque  of  a  coil  of  moment  m 
in  a  magnetic  field  is 

L  -  -mx  B.  (16) 

For  our  geometry  and  parameters  (cf.  Fig.  2), 

L,  =  -mB  cos  6  =  -mB  ,  (17) 

m  m 

if  a  steady  torque  is  assumed.  For  a  flat  coil  with  n  turns  and  area  A  the  moment 
is  given  by 

m  =  nAI,  (18) 

where  I  is  the  current.  The  area  A,  as  before,  is  roughly  the  area  of  the  pole 
pieces  whose  radii  are  r.  Therefore  the  torque  becomes 

L  =  -mrr2IBm  =  -K3I,  (19) 

where  is  defined  as  the  balance -coil  constant.  The  torque  output  is  me¬ 
chanically  coupled  to  the  shaft,  thus  completing  the  feedback  loop.  The  proper 
sign  of  is  obtained  by  correct  polarity  of  the  coil  leads  to  produce  a  negative 
feedback  torque.  The  balance  coil  has  two  taps,  one  at  n  =  10  and  one  at  n  =  100, 
for  control  over  the  system  sensitivity. 

The  Entire  System 

The  block  diagrams  for  the  entire  system  now  become  those  of  Fig.  10. 
Table  1  summarizes  the  various  relationships  thus  far  derived.  The  param¬ 
eters  on  the  right  must  be  found  experimentally  (cf.  next  section):  in  addition, 
several  of  the  constants  (for  example  a,  K^,  and  s^  were  measured  separately 
to  ensure  accuracy. 
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Table  1,  Review  of  parameters. 


J  =  j  m(R2  +  R2) 

=  moment  of  inertia 

a  =  mechanical  damping  term 

m  =  mass  of  shaft 

R^  =  inner  radius  of  shaft 

R  =  outer  radius  of  shaft 

2 

X 

(R  +  rt)  lJ(R  +  rt)2  l 

*1,2  — ziT-i] 1  rp 

=  amplifier  poles 

R  =  4400-ohm  resistor 

r^  =  transformer  resistance 

Lj.  =  transformer  inductance 

C  =  filter  capacitance 

*  2 

K,  =  ntr  B 

3  m 

A  L0 

~  T 

n  =  number  of  turns  on  coil 

r  =  radius  of  pole  pieces  of 

magnetron  magnet 

B  =  field  of  magnetron  magnet 
m 

<l>(s)  _  1 

^e^S^  Js2  +  as  +  K 

.  -  . . 

K  =  shaft  spring  constant 

K2S1S2 

F  Is  +  SjKs  +  sz) 

A  I(s) 

~ETs] 

K.j  =  static  amplifier  gain 
constant 

E(s)  _ 

=  optical  gain  constant 

* 

The  negative  sign  appearing  in  Eq.  19  has  been  put  into  the  feedback  loop 
to  convene  with  standard  notation. 

To  analyze  the  complete  system,  '  the  following  conventional  definitions 
are  made: 


7)  F.  L.  Nixon,  "Principles  of  Automatic  Controls,  "  Prentice -Hall,  New  York, 
1953. 
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(b) 

Fig.  101  Block  diagram:  (a)  general,  (b)  numerical. 


G(s) 

W(s) 


^  (s) 


=  open-loop  transfer  function  (with  no  feedback) 


^  =  closed-loop  transfer  function  (with  feedback) 


The  loop  may  be  physically  severed  by  disconnecting  the  balance -coil  leads  and 
substituting  a  load  resistor  for  the  coil.  The  current  flowing  in  this  load  re  - 
sistor,  multiplied  by  K^,  gives  the  effective  open-loop  output  torque. 

From  the  block  diagram  of  Fig.  10  we  find  that 


G(s)  = 


K1K2K3S1S2 


(Js  +as  +  K)(s  +s^)(s  +s^) 


The  static  gain  of  the  amplifier,  K^,  may  be  varied  at  will  through  the  gain 
control  of  Fig.  4.  The  balance-coil  constant  similarly  is  double -valued 
at  will.  We  assume  operation  on  the  tap  at  n  =  100.  The  remaining  parameters 
of  Eq.  20  are  all  constants.  At  the  sacrifice  of  generality  we  now  substitute 
numerical  values  for  all  parameters  except  K^,  to  permit  a  number  of  simpli¬ 
fications  that  will  ultimately  yield  a  relatively  simple  expression  governing  the 
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Table  2.  Various  numerical  values. 


Quantity 

— 

Symbol 

■ 

Value 

Mass  of  shaft 

m 

3.8xl0-2  kg 

R1 

2. 0x10  2  m 

Outer  radius  of  shaft 

R2 

2.  5  x  10  2  m 

Radius  of  magnetron-magnet 
pole  pieces 

r 

2.  54xl0~2  m 

Magnetron -magnet  field 

Bm 

0.85  weber /m2 

Transformer  inductance 

Lt 

100  hy 

Transformer  resistance 

rt 

7700  ohms 

Total  equivalent  detector 

12, 100  ohms 

resistance 

R 

Optical  gain 

K1 

15.  6  volts /radian 

Static  gain  of  amp 

K2 

variable 

Torque  constant 

K3 

0. 139  n-m/amp 

0.  014  n-m/amp 

Amplifier  poles 

S1 

s. 

20  /sec 

101  /sec 

Moment  of  inertia 

J 

2  x  10  1  kg-m2 

Damping  factor 

a 

2x10  n-m-sec/radian 

Spring  constant 

K 

2x10  ^  n-m /radian 

Number  of  coil  turns 

n 

100  or  10 

behavior  of  the  apparatus  over  the  frequency  range  of  interest. 

Table  2  lists  the  various  numerical  values.  Substituting  these  into  Eq.  20, 

we  get 


G<s)  = 


2.  2x  1010  K, 


(s  +  0.  l)( s  +  10  )(s  +  20)( s  +  101) 


(21) 


We  have  shown  (cf.  p.  5)  that  the  maximum  operating  frequency  is  about  0.4  cps 
(=2.5  radians/sec).  Measurements  showed  responses  up  to  radian  frequencies 
of  about  50.  Hence,  we  may  neglect  the  pole  at  s  =  -10  .  G(s)  may  be  approxi¬ 
mated  by 


2.  2x10  K- 


~  (s  +  0.  l)(s  +  20)(s  +  TUT}"  ' 


(22) 
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Now,  W(s)  may  be  evaluated 

G(s) 


2.2x10  K„ 


W(s)  =  r+T5T5)“  ;2032s ;  2.2xl0k  K.  • 


(23) 


It  is  also  instructive  at  ihi=  point  to  substitute  the  usual  value  for  K?, 

Li 


0. 077  mho. 


" TsT  ii  bj( S T  2!  r  +  j3 8)( s  +  2 .  5  -  j 3 8)  *  (24> 

Further,  if  we  neglect  the  pole  at  s  =  —11 6  (this  is  not  entirely  valid  at  fre¬ 
quencies  around  50  but  useful  nevertheless),  we  get 

vrt  _  1.5x10 

's)  -  T7+‘2:T+''335)(s  +  2.5  -  j38) 


_  k 

(s  +  ad  +  J“0Ks  +  &d  -  jwQ)  * 

where  ad  is  a  damping  term  for  the  response  and  uq  its  natural 
quency. 


(25) 

resonant  fre- 


Experimental  Verification 
Determination  of  Parameters 

The  numerical  values  of  the  parameters  of  Table  1  appear  in  Table  2. 

The  torque  constant  K3  was  evaluated  by  several  different  experimental 
procedures  to  ensure  accuracy,  because  calibration  of  the  instrument  depends 
on  this  value „ 

The  optical  and  static  amplifier  gain  parameters  were  determined  by 
deflecting  the  shaft  by  a  given  angle  and  measuring  E  and  I.  In  this  way  Figs.  11 
and  12  were  obtained. 

The  expression  for  F(s)  of  Eq.  14  becomes 

2020  K2 

(T+' 2 0)(s  +  101)  ’ 


F(s)  = 


(26) 


Relative  output  in  db 


.9 
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of  the  curves  agrees  with  that  predicted  in  that  the  system  is  representable  by 
a  pair  of  complex  conjugate  poles  (cf.  Eq.  25).  The  measured  resonant  fre¬ 
quency  is  about  14,  as  compared  to  38  predicted  by  Eq.  25.  The  measured 
phase  diagram  shows  the  expected  rapid  transition  from  around  0  to  -180  degrees 
at  the  resonance  point. 

To  get  the  experimental  Bode  plots,  a  circular  coil  was  suspended  from 

the  shaft  in  the  field  of  tne  Varian  magnet  and  excited  by  a  sinusoidal  current 

source.  The  presence  of  this  coil  substantially  altered  the  parameters  of  the 

system,  causing  disagreement  with  the  results  of  Eq.  24.  In  particular  the 

coil  increased  the  mechanical  damping  factor  of  electromagnetic  origin  from 
-'1  -3 

2  x10  to  9x10  n-m-sec.  Substituting  this  new  value  of  a  into  Eq.  20  and 
evaluating  the  new  closed-loop  response  function  shows  that  a  resonant  fre¬ 
quency  of  19-  5  radian/sec  and  a  system  damping  factor  of  7.9  are  achieved  for 
the  system.  A  plot  of  the  revised  function  W(s)  in  Fig.  14  shows  good  agree¬ 
ment  with  the  experimental  results.  The  increase  in  the  mass  of  the  shaft  and 
the  equivalent  spring  constant  prove  of  negligible  importance. 


Transient  Analysis 

The  transient  response  of  the  unit  was  measured  by  applying  a  step  of 
torque  to  shaft  (by  brushing  the  shaft  with  a  finger)  and  observing  the  subsequent 
behavior  for  different  gains.  The  system  demonstrated  an  exponentially  damped 
oscillatory  response.  The  step  response  (for  a  unit  step)  may  be  found  from 
Eq.  25.  A  unit  step  has  the  Laplace  transform  l/s.  The  output  in  the  frequency 
domain  is  then 

L  (s)  = 

The  output  response  may  be  inversely  transformed  as 
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Since  «  uq,  Eq.  28  simplifies  to 

LQ(t)  =  — j- (l  -  e  d  cos  GiQt)  =  (1  -  e~2-  5t  cos  38t),  (29) 

~o 

*3 

where  k  =  1.  5x10  (cf.  Eq.  25).  Note  that  for  t  »  l/2.5,  Lq — ►  1,  as  expected. 

•N 

The  step  response  was  plotted  with  a  Sanborn  time -base  recorder.  From 
the  transient  curves  the  following  values  were  found: 


K2 

Q> 

O 

Qd 

0.077 

26.4 

1.53  ' 

0.  135 

31.4 

0.82 

Measured 

0.  147 

34.5 

0.32 

f 

0.  160 

37.7 

0.00 

y 

0.077 

38.0 

2.50 

Calculated 

Again  reasonable  agreement  is  observed  between  theory  and  experiment  for  = 
0.077.  The  values  of  and  are,  of  course,  dependent  on  the  gain  setting. 
For  some  critical  gain  setting,  <  0,  and  the  system  goes  into  an  unstable  con¬ 
dition  of  oscillation.  This  occurred  for  a  gain  of  K?  =  0.160, 

Open-Loop  Considerations 

No  open-loop  data  could  be  taken  because  inherent  open-loop  drift  in  the 
output  stage  of  the  amplifier  would  rapidly  drive  the  system  into  a  nonlinear 
region  of  operation,  as  shown  by  the  fact  that  the  open-loop  system  function 
(cf.  Eq.  22)  has  a  pole  very  close  to  zero  and  hence  a  very  large  open-loop 
gain  for  low  frequencies. 

Stability  Considerations 

An  exaggerated  Nyquist  diagram  of  C(s)/K^  (Eq.  22)  is  shown  in  Fig.  16. 

g\ 

The  Nyquist  criterion  '  states  that  for  stability 

8)  C.  Newton,  L,  A.  Gould,  J.  F.  Kaiser,  "Analytical  Design  of  Linear  Feed¬ 
back  Controls,  "  John  Wiley  and  Son,  New  York,  1957,  p.  302. 
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Fig.  16. 
Nyqaist  plot. 


Z  =  N  +  P  =  0,  (30) 

where  Z  is  the  number  of  right-hand  plane  poles  in  W(s),  N  the  net  number  of 
positive  encirclements  around  a  -l/^C^  point  in  the  Nyquist  polar  plot  of  G(s)/t 
and  P  the  number  of  right-hand  plane  poles  of  G(s).  Clearly  from  Eq.  22,  P  =  0. 
The  Nyquist  diagram  is  broken  into  four  -distinct  regions.  For  region  A  there 
are  no  encirclements  of  the  -l/K^  point;  hence,  N  =  0  and  Z  =  N  +  P  =  0.  Here 
the  system  must  be  stable  where  0  ^  <  0.11  =  l/?.2.  For  region  B,  N  = 

-2  and  thus  the  system  is  unstable  according  to  Eq.  30.  For  region  C,  N  =  -1 
and  again  instability  exists.  For  region  D,  N  =  0;  hence  Z  =  0  and  stability 
is  again  implied.  This  corresponds  to  -l/K^^  1.1x10^.  Actually,  for  any 
<  0  the  system  runs  into  a  highly  nonlinear  region  of  operation  where  the 
coil  form  bangs  against  its  stops;  thus  operation  at  this  value  has  virtually  no 
significance  here. 

Stability  prevails  only  for  values  of  -l/&2  that  lie  entirely  to  the  left  of 
the  polar  plot.  In  the  transient  response  we  found  the  system  to  be  unstable 
for  K2  >  0.16. 


N  is  positive  if  the  plot  is  traversed  in  a  counter-clockwise  direction  for  in¬ 
creasing  values  of  oj. 
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Operation 

Calibration  and  Accuracy 

In  the  torque -measuring  equipment  the  important  parameters  for  cali¬ 
bration  of  the  X-Y  recorder  are  the  balance -coil  constant  and  the  current¬ 
measuring  resistor  Rx-  The  latter  {cf.  Fig.  4)  is  simply  a  small  resistance 
in  series  with  the  balance  coil.  The  voltage  across  R^,  when  plotted  cn  the 
recorder,  gives  the  feedback  current,  which  in  turn  gives  the  feedback  torque 
from  the  relation 

K-E 

h  =  v =  4-r  *  (31> 

X 

where  E^  is  the  voltage  recorded.  may  be  one  of  two  plug-in  resistors. 

There  is  also  a  voltage -dividing  potentiometer  across  Rx  to  permit  a  more 
flexible  adjustment  of  the  recorder  range.  The  angular  circuit  which  monitors 
the  angle  of  the  field  with  respect  to  the  crystal  is  given  in  Fig.  17.  The  Heli- 
pot  is  ganged  directly  to  the  Varian  magnet  and  the  output  voltage  is  directly 
proportional  to  the  angle  of  rotation.  The  circuit  satisfies  the  proper  input 
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requirements  of  the  recorder's  X  axis.  The  two  variable  resistors  permit  ad¬ 
justment  of  the  zero  point  and  compensation  fcr  changes  in  battery  voltage.  The 
accuracy  of  the  angular  circuit  is  limited  to  about  +1%  principally  by  the  read¬ 
ability  of  the  recorder.  Accuracy  for  the  torque  axis  is  conservatively  set  at 
about  2%  for  torques  above  300  ayne-cm.  Below  this  figure,  noise  and  building 

N. 

vibrations  present  increasing  difficulty. 

Errors 

Several  possible  sources  of  error  have  not  yet  been  fully  discussed. 

The  first  is  a  possible  interaction  between  the  rotating  Varian  field  and  the 
torque  equipment.  This  condition  was  found  to  exist  through  an  interaction  of 
the  rotating  field  with  the  phototubes,  but  was  entirely  eliminated  by  mag¬ 
netic  shielding,  which  also  greatly  reduced  noise  pickup. 

A  second  source  of  error  stems  from  a  possible  mis  orientation  of  the 
crystal  plane  on  the  shaft.  X-ray  techniques  are  accurate  to  only  about  +l/2° 
with  our  equipment  in  aligning  a  crystalline  plane,  and  extreme  care  must  be 
used  to  maintain  this  tolerance  in  transferring  the  crystal  from  the  X-ray 
apparatus  to  the  magnetometer.  This  error  is  approximately  proportional  to 
the  square  of  the  angular  deviation  from  orientation,  to  which  must  be  added 
demagnetization-correction  factors  if  the  sample  is  not  ellipsoidal  in  shape. 
Hence,  information  about  the  second-order  anisotropy  constant  may  be  com¬ 
pletely  enshrouded  by  small  errors. 

A  quite  different  problem  arises  if  the  crystal  plane  is  properly  oriented 
but  the  surfaces  of  the  crystal  are  not  parallel  to  this  plane.  Here  a  consider¬ 
ation  of  the  demagnetization  factors  for  the  altered  geometry  becomes  requisite. 
It  would  be  expected  that  at  least  a  2  9 -error  variation  would  be  obtained  for  this 

case. 

* 

A  computation  verified  the  fact  that  the  magnetic  force  exerted  on  the  electrons 
by  the  Varian's  magnetic  field  was  comparable  to  the  electric  force  exerted  by 
the  E  field  of  the  phototube. 
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BEST  AVAILABLE  COPY 


Data  Analysis 


Basic  Method 


Our  essential  approach  to  data  analysis  lies  in  approximating  the  experi- 

9) 

mental  curve  by  a  Fourier  series  '  and  equating  this  to  the  theoretical  torque 
derived  from  the  anisotropy-energy  expansion.  In  this  way  a  considerable 
amount  of  the  information  contained  in  the  curve  is'useth. 

This  procedure  tacitly  assumes  that  there  is  negligible  effect  from  the 
magnetostrictive  anisotropy,  that  the  demagnetization  energy  is  isotropic  in 
nature,  and  that  the  crystal  must  be  saturated  by  a  field  that  is  large  compared 
to  the  demagnetizing  field,  so  that  the  magnetization  vector  lies  in  the  measuring 
plane  and  essentially  parallel  to  the  applied  field.  These  assumptions  are  never 
precisely  attained.  First  the  anisotropy  energy  prevents  the  magnetization  from 
lying  exactly  parallel  to  the  field.  Another  factor  that  may  influence  the  mag¬ 
netization  is,  for  instance,  the  possible  presence  of  an  inhomogeneous  demag¬ 
netizing  field  (compared  to  the  applied  field)  near  the  edges  of  a  disk-shaped 
sample.  This  could  influence  the  torque  so  as  to  make  the  anisotropy  constants 
appear  to  have  a  field  dependence. 

Theoretically,  the  torque  curves  should  be  measured  in  an  infinite  field 
where  the  magnetization  is  directly  parallel  to  the  direction  of  the  applied  field. 
Usually  the  difference  between  high,  practical  fields  and  infinite  fields  is 
negligible;  however,  Kouvel  and  Graham'^  outline  ^  clever  method  that  permits 

extrapolation  to  infinite  fields  when  it  is  not.  This  technique  requires  a  large 

\ 

number  of  torque  curves  to  be  taken  for  increasing  fields.  Fc^r  a  more  complete 
treatment  the  reader  is  referred,  to  the  original  paper. 


9)  R.  F.  Pearson  and  L.  Guilford,  Report  No.  2159,  Mullard  Research  Labo¬ 
ratories,  December,  1957. 

- -  t 

10)  J.  S.  Kouvel  and  C.D.  Graham,  Jr.,  ''Conference  on  Magnetism  and  Magnetic 
Materials,"  AIEE,  Publication  T-91,  February,  1957,  p.  85. 
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Applications  of  the  Instrument 
Anisotropy  Measurements 

Thus  far  measurements  have  been  made  on  several  single  crystals  of 

synthetic  magnetite  (Fe^CO,  nickel,  and  yttrium  iron  garnet  (all  cubic). 

One  magnetite  crystal  v/as  first  measured  in  a  (100)  plane  of  an  ellip- 

-2  3 

soidal  sample  (V  =  1. lx  10  cm  ;  major  axis  about  4  mm).  A  second  was 

-3  3 

ground  into  a  sphere  (Vs  4.4  x  10  cm  )  and  oriented  in  a  (110)  plane.  A  run 
from  liquid -nitrogen  to  room  temperature  was  taken  of  the  latter.  Figure  18 


shows  the  results  along  with  points  for  taken  by  Bickford  '  in  a  resonance 
experiment.  Room-temperature  values  are  summarized  in  Table  3,  The 
values  of  are  generally  consistent  with  those  in  the  literature,  which  run 

C  1 

from  -0.98  to  -1.20x10  ergs /cm  . 

-2  3 

A  room-temperature  measurement  on  a  disk  (V  =  6x10  cm  ;  r  =  9  mm) 

•j-  ^ 

of  nickel  was  run,  giving  values  of  =  -50,  000  and  =  14,  000  ergs/cm  at 

room  temperature.  A  sample -torque  curve  is  shown  in  Fig.  19.  The  value  of 

*  ' 

Grown  in  the  Laboratory  for  Insulation  Research  at  M.I.  T.  by  the  Bridgman- 
Stockbarger  method  (cf.  Tech.  Rep.  49,  December,  1951). 

t 

Purchased  from  the  Virginia  Research  Institute  and  made  available  by  the 
Lincoln  Laboratory  of  M.I.  T. 
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Kj  compares  favorably  with  -49,  00Q 

the  average  of  ten  values  cited  by 

Bozorth.  ^  No  reliable  values  of 

have  been  published. 

Finally  a  temperature  run 

on  an  yttrium  iron  garnet  sphere  * 

-3  3 

(V £=  4.  2 x  10  cm  )  was  made  with 
the  results  shown  in  Fig.  20  along 
with  Dillon1  s*^  data.*  The  anisot¬ 
ropy  differed  in  several  of  the  (110) 
planes,  indicating  that  the  crystal 
was  probably  strained  or  contained 
some  irregular -shaped  impurities. 

A  room-temperature  measure ¬ 


Fig.  19.  Torque  vs.  angle  curve  for  nickel 
in  the  (110)  plane  at  20°C. 


ment  on  a  single -crystal  disk  with 

a  greater  volume  (V  =  2.8x10  ^  cm3)  gave  reasonable  agreement  with  the  previous 
results. 

Table  3.  Room-temperature  anisotropy  constants  (in  ergs/cm  ). 


K1 

K2 

Magnetite  ellipsoid 

-1.13  xlO5 

- 

Magnetite  sphere 

-1. 12  x  105 

(6  +  3)  x  104 

Nickel  disk 

-5. 0  x 104 

(1.4  +  0.  5)  x  104 

Yttrium  iron  garnet  sphere 

-6. 1 x  103 

(1.0  +  0.  5)  x  104 

Yttrium,  iron  garnet  disk 

-6.  5xl03 

(0.5  +0.5)xl04 

* 

Grown  in  a  PbO-PbF^  flux  by  a  process  similar  to  that  described  by  Nielsen, 
Suppl.  to  J.  Appl.  Phys.  3R  No.  5,  51S  (I960).  The  sample  was  obtained  from 
Air  Force  Cambridge  Research  Center. 

11)  J.  F.  Dillon,  Phys.  Rev.  105,  759  (1957). 


To  find  K]^  from  Dillon’s  data  given  in  terms  of  K^/Mg,  we  used  values  of 
found  by  Pauthenet,  J.  Appl.  Phys.  29,  253  (1958). 
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Rotational  Hysteresis 

Since  the  torque  magneto¬ 
meter  measures  minute  torques, 
it  can  also  be  used  to  3tudy  rota¬ 
tional  hysteresis.  If  a  sample  is 
suspended  in  a  uniform  magnetic 
field  and  rotated  in  this  field  (or 
the  field  is  rotated  around  it),  the 
rotational  hysteresis  is  defined 


Fig.  20.  Kj  vs.  temperature  for 
'  -  yttrium  iron  garnet. 


as  the  energy  loss  per  radian  per  unit  volume: 

/.2tt  \ 

Wr  =  -27P  f  L(°)d9'  (32) 

•'o 

where  V  is  the  volume  of  the  sample  and  the  torque  is  given  in  dyne-cm.  The 
rotational  hysteresis  energy  loss  represents  that  amount  of  work  that  must  be 
done  on  the  crystal  in  order  to  rotate  the  magnetization  through  1  radian. 

It  usually  appears  in  the  crystal  as  thermal  energy. 

The  rotational  hysteresis  of  a  polycrystallinc  yttrium  iron  disk  sample 
as  a  function  of  the  measuring  field  is  shown  in  Fig.  21. 


Induced  Anisotropy 

In  certain  polycrystalline  materials, 
such  as  yttrium  iron  or  lutetium  iron  garnet, 
a  uniaxial  anisotropy  may  be  annealed  into 
the  specimen  by  cooling  in  a  large  field. 

For  polycrystnlline  YIG  n  uniaxial  anisot¬ 
ropy  of  about  900  ergs/cm^  was  induced  by 
annealing  in  a  6-kg  field  from  room  to  liquid- 
nitrogen  temperature.  This  represents 
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Fig.  21.  Rotational  hysteresis 
of  polycrystalline 
yttrium  iron  disk. 
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Fig.  22.  Field -induced  anisotropy  in  polycrystalline  yttrium  iron  garnet. 

roughly  5%  of  the  single -crystal  value.  For  I_1G  a  value  of  about  3500  ergs /cm 
(or  50% of  the  single -crystal  value)  was  similarly  obtained.  By  altering  the  di¬ 
rection  of  the  applied  field  when  at  77°K,  the  axis  of  symmetry  can  be  rotated 
(Fig.  22). 

Decay  Mechanism 

It  has  also  been  found  that  a  relaxation  process  is  associated  with 
materials  exhibiting  the  previous  anneal  effects  that  can  be  observed  with  the 
torque  magnetometer.  This  may  be  seen  by  annealing  the  sample  as  before 
and  subsequently  altering  the  direction  of  the  applied  field  by  45  degrees  to 
coincide  with  the  position  of  maximum  torque  as  dictated  by  the  uniaxial  ani¬ 
sotropy.  There  will  be  a  torque  due  to  the  "initial”  anisotropy,  which  will 
decay  with  time  approximately  as  an  exponential  as  the  easy  direction  of  mag¬ 
netization  relaxes  over  to  the  new  position.  This  relaxation  process  for  YIG 
has  a  time  constant  of  about  30  seconds  (Fig.  23). 
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Fig.  23. 

Torque  relaxation  curve  for 
yttrium  iron  garnet. 


Conclusion 

The  torque  magnetometer  has  been  developed  to  the  point  where  it  may 
be  used  for  routine  measurements  of  magnetic  anisotropy,  rotational  hysteresis, 
and  anneal  effects  over  a  wide  temperature  range,  to  be  ultimately  extended  to 
liquid  helium. 

Acknowledgment 

The  author  is  indebted  to  the  many  members  of  the  Laboratory  for 
Insulation  Research  who  offered  their  skillful  assistance  and  advice  through¬ 
out  all  phases  of  this  project.  In  particular  he  would  like  to  mention  Profes¬ 
sor  D.  J.  Epstein  and  Mr.  B.  Frackiewicz  for  their  manyihelpful  suggestions 
and  discussions. 


